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ABSTRACT 



We present an analysis of the Internal Shock Model of GRBs, where gamma-rays 
■ are produced by internal shocks within a relativistic wind. We show that observed 

. GRB characteristics impose stringent constraints on wind and source parameters. 

^ i We find that a significant fraction, of order 20%, of the wind kinetic energy can be 

converted to radiation, provided the distribution of Lorentz factors within the wind has 
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a large variance and provided the minimum Lorentz factor is > T± « lO'^'^L^^, where 



L = 10^^L52ergs ^ is the wind luminosity. For a high, > 10%, efficiency wind, spectral 
CN ■ energy breaks in the 0.1 to 1 MeV range are obtained for sources with dynamical time 

^ ■ R/c ^ 1 ms, suggesting a possible explanation for the observed clustering of spectral 

. break energies in this range. The lower limit Tj- to wind Lorenz factor and the upper 

limit ~ l(iZ/10^cm)~^/^ MeV to observed break energies are set by Thomson optical 
depth due to pairs produced by synchrotron photons. Natural consequences of the 
model are absence of bursts with peak emission energy significantly exceeding 1 MeV, 
and existence of low luminosity bursts with low, 1 keV to 10 keV, break energies. 
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iJ ■ 1. Introduction 

. P.. 

The widely accepted interpretation of the phenomenology of 7-ray bursts (GRBs) is that the 
observable effects are due to the dissipation of the kinetic energy of a relativistically expanding 
wind, a "fireball" (see ^ eszaros 1995| and piran 19"96| for reviews). The discovery of GRB afterglow 



emission in X-ray (Costa et al. 1997), optical (van Paradijs et al. 1997) and radio wavelengths 
(Frail et al. 1997) confirmed (Waxman 1997, Wijers, Rees & Meszaros 1997), standard model 
predictions of afterglow (Paczyhski & Rhoads 1993, Kats 1994, Meszaros & Rees 1997, Vietri 
1997) that results from the collision of the expanding fireball with surrounding medium. 

Both the spectrum and temporal dependence of afterglow emission are consistent with 
synchrotron emission of electrons accelerated to high energy at the shock wave driven by the 



^Osservatorio di Arcetri L. E. Fermi 2 50100 Firenze, Italy 

^Department of Condensed Matter Physics, Weizmann Institute, Rehovot 76100, Israel 



- 2 - 



fireball into its surrounding medium. The rapid temporal variability of 7-ray emission, i.e. of 
the GRB itself, imply, on the other hand, that the GRB is produced by internal shocks within 
the expanding wind ( Woods &: Loeb 1995| , [Sari fc Piran 1997 ; Daigne & Mochkovitch, 1998; 
Panaitescu, Spada & Meszaros 1999). Both synchrotron and inverse-Compton emission from shock 
accelerated electrons have been proposed as the GRB emission mechanism. However, synchrotron 
emission is favored if the fireball is required to be "radiatively efficient," i.e. if a significant 
fraction of the fireball energy is required to be converted to 7-rays (e.g. perishev et al. 20001 ). 



Observations of GRB970508 afterglow implies a high radiative efficiency for this GRB ( Frail 



Waxman &: Kulkarni 2000|) , and X-ray afterglow observations imply that high radiative efficiency 
is required for most GRBs ( [Freedman fc Waxman 2000 ). 



Observations therefore suggest that the GRB emission is due to synchrotron emission of 
electrons accelerated by internal fireball shocks. This model faces, however, three major difficulties. 
First, it is generally argued that if 7-ray emission is due to internal shocks within the fireball, then 
only a small fraction, < 10^^, of fireball energy is converted to 7-rays (Daigne & Mochkovitch, 
1998; [Kumar 1999| , [Kumar fc Piran 1999| , jSpada, Panaitescu fc Meszaros 2000[) . This is in conflict 
with the high radiative efficiency implied by afterglow observations. Second, the 7-ray emission of 
more than 60% of the GRBs observed by BATSE detectors is peaked at photon energies in the 
range of 50 keV to 300 keV (Brainerd et al. 1999). It is generally claimed that clustering of peak 
emission energy in this narrow energy range requires fine tuning of fireball model parameters, 
which should naturally produce a much wider range of peak emission energies. Finally, there is 
evidence that, at least in some cases, the GRB spectra at low energy is steeper than expected for 



synchrotron emission (Preeece et al. 1995, Frontera et al. 200C). 



The steep low energy spectra has led several authors to consider inverse-Compton emission 
as the source of 7 radiation (e.g. Lazzati, et al 2000). However, steep low energy spectra may 
also be accounted for within the context of the synchrotron emission hypothesis by the effects of 
inverse-Compton suppression at low energies ( Perishev et al. 2000] ) and by the contribution to 
observed 7-ray radiation of photospheric fireball emission (e.g. Meszaros & Rees 2000). As we 
show here (see discussion in §4), the contribution of latter effect is expected to be significant. 

The main goal of the present paper is to determine the constraints imposed by GRB 
observations on the relativistic fireball wind parameters, and hence on the underlying source 
producing the wind, under the assumption that the observed radiation is due to synchrotron 
emission produced in internal wind shocks. In particular, we address the questions of whether, 
and under what conditions, high radiative efficiency can be obtained, and whether the clustering 
of peak emission energies can be naturally explained by the model. 



It has been demonstrated in earlier work (e.g. Kobayashi et al. 1997; Beloborodov 2000), 
that high radiative efficiency may be obtained in the internal shock model provided the variance 
of the wind Lorentz factors distribution is large. Kobayashi et al. (1997) found ~ 10% efficiency 
for uniform Lorentz factor distribution with maximum to minimum Lorentz factor ratio ~ 10, 
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and Beloborodov (2000) demonstrated that the efficiency may approach 100% for non-uniform 
distributions with wide range of Lorentz factors^. In the analysis presented here we consider 
several issues which have not been addressed in earlier work: The constraints imposed on the 
model by the observed peak emission energy distribution, the effect of optical depth due to 
pairs produced within the fireball wind, and the upper limit imposed on the maximum wind 
Lorentz factor by the acceleration process. The analysis presented here is also more general, as we 
study the dependence of radiative efficiency on a wider set of models parameters. We find that 
the radiative efficiency is limited to values well below the 80%-90% derived by Beloborodov. The 
inclusion of pair optical depth effects is essential in determining both the efficiency and the 
peak emission energy. 

High radiative efficiency clearly requires that a large fraction of the internal energy generated 
by internal shocks be carried by shock accelerated electrons, i.e. that the electron and proton 
energy densities be close to equipartition. Moreover, high radiative efficiency requires the magnetic 
field energy density to be not far below equipartition as well (e.g. Derishev et al. 2000). Thus, 



we assume in most our calculations that electrons, protons, and magnetic field energy densities 
are close to equipartition, and focus on the study of the dependence of observed 7-ray flux and 
spectrum on wind luminosity and on wind dynamical parameters (source size, variability time, 
Lorenz factor and mass distribution of shells within the wind). For completeness, we demonstarte 
using our calculations (in §3.2) that near equipartition is required in order for model predictions 
to be consistent with the observed peak emission energy distribution. 

The model outline is presented in §2. The results of a numerical analysis of the model are 
presented and discussed in §3. The main conclusions are summarized in §4. 



2. Outline of the model 

In the fireball model of GRBs, a compact source, of linear scale Rq ~ 10^ cm, produces a 
wind characterized by an average luminosity ~ lO^^ergs^"*^ and mass loss rate Af = L^/r]c^ 
{Rq 10^ cm corresponds to three times the Schwarzschild radius of a non-rotating, solar mass 
black hole, and ~ lO^^ergs"^ is the characteristic luminosity inferred from observations). At 
small radius, the wind bulk Lorentz factor, F, grows linearly with radius, until most of the wind 
energy is converted to kinetic energy and F saturates at F ~ 7/ ~ 300. r] ^ 300 is required to reduce 
the wind pair-production optical depth for observed high energy, > 100 MeV, photons to less than 
unity. Ifi] > 7]^ ^ {arLyj / ATTiTipC' RqY/"^ = 2 x 103(L^,52/i?o,6)^''^, where = lO^^^^^^^erg s"^ 
and i?o = 10^i?o,6 cm, the wind becomes optically thin at F ~ ry^, < r;, and hence acceleration 
saturates at F ~ ry,,, and the remaining wind internal energy escapes as thermal radiation at 



^Kobayashi & Sari (2001) similarly find ~ 60% efficiency for non-uniform distribution of Lorentz factors over the 
range 10 to lO"*. 
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r-^ 1 MeV temperature. Variability of the source on time scale ty, resulting in fluctuations in the 
wind bulk Lorentz factor T on similar time scale, then leads to internal shocks in the expanding 
fireball at a radius 

Ri r^ct^ = 3 X lO^^r^t^ _3 cm, (1) 

where F = lO^T^;, ty = 10~'^t^^_3 s. If the Lorentz factor variability within the wind is significant, 
internal shocks reconvert a substantial part of the kinetic energy to internal energy. It is assumed 
that this energy is then radiated as 7-rays by synchrotron (and inverse-Compton) emission of 
shock-accelerated electrons. 

In this work, we use an approximate model of the unsteady wind described in the preceding 



paragraph [following Kobayashi et al. 1997 ; Panaitescu, Spada, &: Meszaros 1999; Spada, 
Panaitescu & Meszaros 2000 (SPMOO)]. The wind evolution is followed starting at radii larger than 
the saturation radius, i.e. after the shells have already reached their final Lorenz factor following 
the acceleration phase, and the GRB photon flux and spectrum resulting from a series of internal 
shocks that occur within the wind at larger radii are calculated. We consider both synchrotron 
and the inverse-Compton emission, and take into account the effect of pair production. 

We approximate the wind flow as a set of discrete shells. Each shell is characterized by four 
parameters: ejection time tj, where the subscript j denotes the j-th shell, Lorentz factor Tj, mass 
Mj, and width A^. Since the wind duration, tw ~ 10 s, is much larger than the dynamical time 
of the source, td ~ Ro/c, variability of the wind on a wide range of time scales, td < ty < t^, 
is possible. For simplicity, we consider a case where the wind variability is characterized by a 
single time scale > ty > td, in addition to the dynamical time scale of the source td and to the 
wind duration tyj. Thus, we consider shells of initial thickness Aj = ctd = Rq, ejected from the 
source at an average rate t~^. The shell thickness may increase with time, due to variations in the 
expansion Lorenz factor across the shell. The Lorenz factor distribution across the shell depends 
on the entropy distribution, which is determined by the details of the ejection from the source. 
We therefore consider two extreme cases: uniform Lorenz factor across the shell, in which case the 
shell thickness is time independent, and order unity variation of LF across the shell, AF/F ~ 1, in 
which case the shell thickness is given by Aj = max{Ro,R/Tj). 

We assume that the Lorenz factor (LF) of a given shell is independent of those of preceding 
shells, and consider 4 different LF distributions: uniform, modulated, bimodal, and lognormal. In 
the uniform distribution case, the shells' LFs are randomly drawn from a uniform distribution over 
the range F^ to Tm, where F^ < F^/ are time independent. In the modulated case, Tm varies 
as sm'^{2TTtj/ty,). In the bimodal case, the LFs are drawn from a bimodal distribution, F-,- = F^ 
or Tj = Tm with equal probability. Finally in the lognormal case, shell LFs are drawn from a 
lognormal distribution with an average (F) = 1000 and {T'^y/'^/{T) = 4, truncated at F > ry^. 

The time intervals tj+i — tj are drawn randomly from a uniform distribution with an average 
value of ty. The shell masses Mj are chosen to be either independent of j, or inversely proportional 
to Tj. That is, we consider the two qualitatively different possibilities: equal shell mass and 
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equal shell energy. The total mass is determined by the constraint J2jLi ^j^jC^ = Lyjtw, where 



Once shell parameters are determined, we calculate the radii where collisions occur and 
determine the emission from each collision. We assume that following each collision the two 
colliding shells merge and continue to expand as a single shell. Under this assumption, the internal 
energy carried by the merged shell following the collision is 



is the expansion LF of the merged shell and Mi 2, Ti 2 are the masses and the LFs of the colliding 
shells. The dynamical efficiency e^, defined as the fraction of total kinetic energy converted to 
internal energy, increases with the ratio of the fast and slow shell LFs, and, for a given ratio, is 
maximized for equal masses. The first collisions, for which the differences in the shell LFs are 
largest, are the most efficient. 

In each collision a forward (FS) and a reverse (RS) shock are formed, propagating into 
forward and backward shells respectively. The plasma parameters behind each shock are 
determined by the Taub adiabat, requiring continuous energy density and velocity across the 
contact discontinuity separating the two shells (Panaitescu k. Meszaros 1999). Since the energy 
density in the downstream regions of both shocks is similar, we divide the internal energy between 
the reverse and forward shock according to the ratio Ef/E^ = Af/kAr, where A^j are the shell 
thicknesses prior to the collision and A: is a factor of the order of few that takes into account the 
compression ratios of the two colliding shells. 

The energy released in each shock is distributed among electrons, magnetic field and protons 
with fractions eg, cb and 1 — (eg + es) respectively. We assume that electrons are accelerated by the 
shocks to a a power-law distribution, dne/dje oc 7"^ for particle Lorentz factors 7inin < 7e < 7max- 
7max is the electron Lorenz factor at which the synchrotron cooling time equals the acceleration 
time, estimated as the electron Larmor radius divided by c, while 7inin is determined by the 
requirement that the energy carried by electrons be equal to teEin- 

The GRB energy spectrum and flux is derived by summing the contributions of individual shell 
collisions. For each collision, synchrotron and inverse-Compton emission by the shock-accclcrated 
electrons is calculated. In order to achieve high radiative efficiency, electrons must radiatively 
loose most of their energy on a time scale shorter than shock shell crossing time. Under this 
condition, the synchrotron spectrum may be approximated as 



Ein = (MiFi + MaFs - (Mi + M2)^)c^ 



(2) 



where 




(3) 




(4) 
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where Vgy is the characteristic synchrotron emission frequency of electrons with 7 ~ 7min- The 

inverse Compton spectrum has a similar shape but is shifted to higher energy as described 
in SPMOO. In determining the inverse-Compton flux, we take into account the Klein-Nishina 
suppression of the inverse-Compton cross section. 

The photons radiated can be scattered by the electrons within the shell. The Thomson optical 
depth tt for such scattering is evaluated by taking into account the electrons that were accelerated 
but have cooled radiatively while the shock crossed the shell, and those within the yet un-shocked 
part of the shell. The optical depth to Thomson scattering may be increased significantly beyond 
the value derived taking into account cooled shell electrons, due to the production of pairs. 
This contribution has not been taken into account in previous analysis of the internal shock model, 
since for a uniform distribution of Lorentz factors, and under the hypothesis of equipartition, 
the break energy, hiysy, in the shell co- moving frame is well below the pair production threshold. 
However, since the photon spectrum extends to high energy as a power law with spectral index 
—p/2 Ri —1, there is an equal number of photons per logarithmic energy intervals, and there may 
exist a large number of photons beyond the pair production threshold. The pairs produced by 
these photons may contribute significantly to the Thomson optical depth. 

In order to take the effect of pair production into account, we determine for each collision the 
photon energy epairs for which the pair production optical depth T^,y equals unity. All the photons 
that have energies higher than max(meC^, Cpairs) (measured in the shell frame) are assumed to 
form pairs, leading to a suppression of the high energy photon tail. We assume that the energy of 
photons that undergo pair production is converted to sub-relativistic pairs, which are taken into 
account in determining the Thomson optical depth. 

We assume that a fraction (1 — e^'^'">^)/rr,/, where r^j are the optical depths of the RS and 
FS, of the photons produced by each shock escapes and reaches the observer (The RS flux is also 
absorbed by the forward shell, and is thus further reduced by a factor e~'^f). The energy of the 
absorbed photons, as well as the internal energy which has not been converted to radiation, are 
converted back to ordered kinetic energy by adiabatic losses during the shell expansion. During 
this process, we take into account the increase in shell thickness by assuming this thickness grows 
as the shell's (comoving frame) speed of sound. 

Our assumption, that following a collision two shells merger into one shell which expands 
with a single (uniform) Lorenz factor, leads to an overestimate of the reduction due to collisions 
of shell Lorenz factor variance. Clearly, the expansion of the shocked shells due to the pressure 
produced by the shocks results in the leading edge of the two shells moving faster (in the observer 
frame) than the trailing edge. It is therefore clear that the variance in Lorenz factor distribution 
following the collision is larger than obtained under our assumption of a "complete" merger. This, 
in turn, implies that a more detailed calculation of the shell merger process, which may require 
dividing the merged shells into several distinct shells following the collision, will enhance the 
overall efficiency of kinetic to thermal energy conversion. The magnitude of such enhancement 
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depends on the density and Lorenz factor distributions within the shells prior to the collision, 
and is expected to be of order unity. We therefore do not consider this effect in the approximate 
analysis presented in this paper. 

3. Results and discussion 

In this section we determine the dependence of the wind radiative efficiency e^, defined as 
the fraction of wind energy converted to radiation, and peak emission energy e^, the photon 
energy at which the maximum of vfi, is obtained, on wind model parameters. As explained in 
the introduction, we adopt electron and magnetic field energy fractions close to equipartition and 
analyze the dependence of and Sp on wind luminosity, variability time, source size, and on 
the distribution of wind Lorentz factors and shell masses. Wc use ee=0.45, cb = 0.1, and p = 2 
throughout the paper, except in the calculations presented in Fig. 4, where smaller eg and 
values are assumed in order to demonstarte that near equipartitrion is required by observations. 
We first discuss in §3.1 the various Lorenz factor distributions described in §2. We then present a 
detailed discussion of the bimodal case, which best reproduces observed GRB characteristics, in 
§3.2. 

3.1. Comparison of various Lorenz factor distributions 

Fig. 1 shows the results of numerical simulations of the model described in the previous 
section. Each panel shows versus Sp for one of the LF distributions described in §2. Different 
points within each panel correspond to different choices of Tm, Tm and t^. Results are shown for 
10 < Tm < Fm < 2500 and 10~^ s < < 1 s. We have assumed equal mass shells, Rq = 10^ cm, 
= lO^^ergs"^ and source redshift 2; = 1 for all calculations. Results are presented for the two 
extreme assumptions on shell width evolution with radius described in §2, i.e. for both constant 
width A = Rq and maximal expansion, A = max(i?o, R/V^). 

Based on Fig.l, a uniform LF distribution can be ruled out since the radiative efficiency is 
small, < 5%. The low radiative efficiency in this case is due to low dynamical efficiency, i.e. to 
the fact that only a small fraction of the kinetic energy is converted to internal energy in collisions. 
In order to increase the dynamical efficiency it is necessary to cither enhance the variance of 
the LF distribution of colliding shells at the first set of collisions (as is the case for bimodal or 
lognormal distributions) or to keep the variance moderate over a wind range of collision radii (as 
is the case for the modulated case). 

In the modulated case the maximum dynamical efficiency is 40%. The first collisions remove 
the initial random differences and the merged shells have LFs near the ejection average value 
(F) = (F^ + Fm)/2. If the wind is homogeneous (F) is similar for all shells, resulting in a steady 
decrease of during the wind expansion. If the range of LFs is varying on time scale of order 
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of tyi], (r) reflects the initial modulation of Tm and larger radii collisions that are dynamically 
efficient are still possible. Thus, a modulation in the wind can solve the efficiency problem 
(e^ ~ 10%). However, the peak emission energy in this case is well below the BATSE range, since 
most of the emission originates at large radii, much larger then the radius of the photosphere R±, 
the radius at which the wind becomes optically thin to Thomson scattering by pairs [see the 
discussion preceding Eq. (||) in §3.2]. 

The dynamical efficiency in the case of a bimodal LF distribution can reach 50%, since 
the Lorentz factors of the colliding shells are very different. Contrary to the case of modulated 
ejection, the collisions at the photospheric radius R± dominate wind emission, and a region of the 
parameter space exists, where model radiative efficiency and peak emission energy are consistent 
with observations, e.y > 10% and Ep ~ 100 keV. 

The case of a lognormal LF distribution can be considered as an intermediate case between 
the random and bimodal ones. In this case the maximal radiative efficiency expected is of the 
order 10-15%, much lower than the 80-90% value derived in Beloborodov (2000). This is due to 
the pair optical depth, which sets a lower limit of F^ ^ 100 for an optically thin wind (see §3.2 
below), and to the constraint imposed on the maximum LF by finite source size, Fa./ < t]* ~ lO'^'^. 
These constraints limit the variance of the LF distribution, and hence the dynamical efficiency. 

Considering shells of equal energy, rather than of equal mass, the main results do not change. 
Uniform LF distribution is still characterized by low efficiency and modulation leads to low peak 
emission energy. In the case of a bimodal distribution, the efficiency is lower for shells of equal 
energies than for shells of equal mass, due to an increase in wind optical depth. The LF of the 
merged shell, given by Eq.(^), is reduced in the equal energy case by a factor ^/T^^/Tm with 
respect to the equal mass case. Thus, the average wind LF is lower, implying smaller collision 
radii and higher opacity. 

^From the analysis of Fig. 1 we conclude that a bimodal distribution of shell LFs, that is, 
LF distribution with large variance, and equal shell mass are favored by observations. Note, 
that this conclusion is independent of the assumption of near equipartition, using which the 
above calculations were performed, since for electron and magnetic field energy fractions below 
equipartition the peak emission energy would be lower than obtained above [see also §3.2, Fig. 
4 and Eq. (^), below]. In what follows we study the bimodal LF case in more detail, in order to 
examine the dependence of observable characteristics on wind parameters. 

3.2. Bimodal LF distribution 

Figures 2 and 3 present the dependence on F^ and on of peak emission energy, Sp, and 
radiative efficiency, e^, respectively, for Tm = 2500, = lO^^ergs"^, Rq = 10^ cm and equal shell 
masses. Results are presented in Fig. 2 for the two extreme assumptions on shell width evolution 
with radius, A = Rq and A = max{RQ, R/T'^). The efficiency contour plot, shown in Fig. 3, is 
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similar for both cases. 

The dependence of Sp on parameters can be understood considering the dependence of the 
dynamical efficiency and of the optical depth on Tm and ty. For high values of Tm (> 300), both 
reverse and forward shocks are optically thin. Decreasing Tm leads to a decrease in the initial 
radius of the collisions, increasing the efficiency up to 15% and the peak emission energy up to 
0.1-1 MeV. At lower values of r^, the efficiency decreases and steep breaks are observed in the Sp 
contour plot. These breaks can be understood considering the difference between the FS and RS 
efficiencies, due to different comoving densities of the two colliding shells {Ef/Er ~ Pc,r/Pcj)- In 
the case of maximal shell expansion, A = R/T"^, the ratio Ef / E^ ~ Fr/Fj, and the FS dominates 
the emission. Thus, the break in the ep contour plot (left panel of Fig. 2) corresponds to the 
region of F^ and values where the FS, and consequently also the RS, become optically thick 
due to pair production (F^ = 50 to F^ = 130 for ty varying between 1 s and 10~^ s respectively). 
A small decrease in F^ results in collisions below the photosphere, leading to a steep decrease 
in peak emission energy, to a value corresponding to the second set of collisions, which occur at 
larger radii. These collisions dominate the emission and are characterized by a lower efficiency 
and a lower peak emission energy. In the case of a constant shell width A = Rq, the ratio 
Ef/Er ~ Tj/rr) and the RS dominates the emission. Thus, there are two breaks in the peak 
energy contour plot (right panel of fig 2). The first corresponds to the values of F^ and ty where 
the RS becomes optically thick due to pair production and the FS starts to dominate the emission. 
This break represents the shift of the peak energy from the RS value to the FS one. The other 
occurs in the parameter region where the FS becomes optically thick, and corresponds to the shift 
of the peak energy from the first to the second set of collisions. 

The contour plot of the radiative efficiency in Fig. 3 shows that values of higher than 10% 
and peak emission energy between 0.1 and IMeV are reached at similar regions in the F^- ty 
plane, in accordance with the correlation between the peak emission energy and radiative efficiency 
shown in Fig. 1. High efficiency, > 10%, is obtained only for large Tm, Tm ~ r/^, « 2 x 10^, i.e. 
for values close to the maximum allowed by the finite source size. For lower values of Tm, the 
dynamical efficiency decreases leading to lower Sp and e-y. 

In Fig. 4 we present results of calculations similar to those presented in Fig. 2, for eg and €b 
values below equipartition. The results demonstarte that the electron and magnetic field energy 
fractions can not be far below equipartition, eg ^ 10^^ and 10~^ are required in order to 

obtain peak emission energies consistent with observations. These results are consisitent with our 
analytic estimates, Eqs. (|8|) and (P) below. 

The contour plot shown in Fig. 2 demonstrates that values of £p larger than ~ 1 MeV can 
not be obtained. This result can be understood using the following arguments. For a collision of 
shells of thickness A at radius R, we have 




(5) 
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The maximal value of Ep is obtained for collisions at the smallest radius R for which the wind is 
optically thin. The radius of the Thomson photosphere due to the electrons present in the original 
fireball is given by 

axM II 1/2 ,1/2 1^-1/2 

^ V 4^ ™ '"'52 K-z > (6) 

where M is the shell mass, and we have approximated M = Lyjtv/(?{Xm + rM)- -^±5 the Thomson 
photosphere radius due to pairs resulting from pair production interaction of synchrotron 
photons, can be estimated by assuming that a significant fraction, ~ 1/2, of the radiative energy 
is converted to pairs, as typically is the case. The number density of is given in this case by 
n± 7=i eeEiri/iSTTmeC^R'^T^A) and 



R± 



y SvrmeC^ V ^ 

inl3^^ ^1/2 rl/2 ,1/2 /-p -p \-l/4 
lU cm -L^ 52 _3 (1 M,3i m,2j ' • 



Comparing R± and Rt, we conclude that the optical depth is typically dominated by the pairs. 

For fixed and Tm, the maximal peak energy is obtained for = F^i for which the 
collision radius Ri ^m'^t^ equals the photospheric radius R±. Solving for Tm±, 

F ^ ~ 5 X f2/9 p-1/9 .-2/9 .2/9 

and substituting in Eq. (^) we find 

^max _ . V ^1/2 4/3 y.4/3 . -1/2 1/6 .-1/6 

~ 2MeV e^ eg' i m,3 ^6 -3 ^«),52 > l^J 



where we have used 7min = \/(rM/rm)eemp/melog(7max/7mm) with 7max = 1007min- The 

dependence of on wind luminosity is demonstrated in figures 5 and 6. Under the assumption 
of fixed shell width, A = Rq, e™^^ is weakly dependent on L, as indicated by Eq. (^. If the 
maximum Lorenz factor scales as F^ r]^ oc L^/^, rather than being independent of L, the slow 
decrease with L of e^^^ at fixed t^, e^^^ oc L~^/^, is replaced by a slow increase with L 

« 4MeV ef Cs L%. (9) 

The dependence on luminosity becomes stronger under the assumption of maximal expansion. 
Substituting A = R/V'^ in Eq. (D, we find e™^^ oc L ^1^^ . The stronger dependence is also 
apparent in figures 5 and 6. However, allowing for a scaling Tm ^ rj^, (x L^/^, the peak emission 
energy dependence on L becomes e™'^^ oc L^/^^ ^ similar to the dependence shown in Eq. (^). 



4. Conclusions 

We have analyzed a model of GRBs, in which a compact source of linear scale Rq produces 
a wind characterized by an average luminosity L.^ and mass loss rate M = Lwjr](? . At small 
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radius, the wind bulk Lorentz factor, T, grows linearly with radius, until most of the wind 
energy is converted to kinetic energy and T saturates at F ~ 77. Variability of the source 
results in fluctuations in the wind saturation Lorentz factor F, leading to internal shocks in the 
expanding wind. These shocks reconvert a fraction of the kinetic energy back to internal energy, 
which is assumed to be radiated as 7-rays by synchrotron (and inverse-Compton) emission of 
shock-accelerated electrons. Since the wind duration, ~ 10 s, is much larger than the dynamical 
time of the source, td ^ Ro/c, variability of the wind on a wide range of time scales, < < tw, 
is possible. For simplicity, we have assumed that in addition to t^, which determines the initial 
shell thickness A = ctd ^ Rq, and t^, the wind is characterized by a single time scale > t^, 
which determines the shell ejection rate t~^. We have addressed the questions of whether, and 
under what conditions, high radiative efficiency consistent with observations can be obtained, 
and whether the observed clustering of peak emission energies can be naturally explained by the 
model. 

We have shown that a significant fraction, ~ 15%, of the fireball energy can be converted 
to radiation. As pointed out at the end of §2, our simplified treatment of post-collision shell 
evolution, assuming complete merger of shells, leads to an overestimate of the reduction of the 
variance of colliding shells' Lorenz factors with wind evolution. A more detailed calculation of 
the merger process will lead to an enhancement of the radiative efficiency by a factor of order 
unity. The exact value of the enhancement factor will depend, however, on the unknown internal 
shell structure. Adopting simplifying assumptions regarding the shell structure and post-collision 
evolution, Kobayashi & Sari (2001) have recently pointed out that including post-collsion evolution 
effects may lead to radiative efficiency exceeding Ce (note, that this may be the case also when 
post-collision evolution is neglected). Although this effect may relax somewhat the constraint 
imposed on eg by the requirement of high efficiency, we have shown here (see following paragraph), 
that ee is constrained to be of order unity by the observed 7-ray spcetra. 

In order to obtain high radiative efficiency and peak emission energy ~ 1 MeV, the minimum 
radius Ri at which internal collisions in the expanding wind occur is required to be similar to 
R±, the radius where the Thomson optical depth due to pairs produced by shock synchrotron 
emission equals unity [see Eq. (0)], and a large variance (compared to the mean) of the colliding 
shells' Lorenz factor (LF) distribution is required. In addition, the electron and magnetic field 
energy fractions should be close to equipartition, with eg ^ 0.1 and es ^ 0.01, in order for 
model predictions to be consistent with the observed peak emission energy distribution [see 
Fig. 4 and Eq.(^)]. The constraint Ri '-^ R± is equivalent to a constraint on the minimum LF, 
Fm, of expanding shells, F,^ ~ F^-i- 

~ 102-5(L^,52/t,,_3)2/9, where = W^^L 52ergs and 
= 10~^tt,^_3 s [see Eq. (^]. Large variance in the LF distribution of colliding shells than requires 
a non-uniform LF distribution (e.g. truncated log-normal or bimodal distributions, see Fig. 1) 
with Tm, the maximum LF of wind shells, close to the upper limit set by the shell acceleration 
process, T < rj^ ^ [gtL^ / ATrnipC^ RqY/^ = 2 x 10^(L^„,52/-Ro,6)^^^, where Rq = 10^i?o,6 cm. 

Our results do not agree with the very high, ~ 80%, radiative efficiency values, obtained in 
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the calculations of Beloborodov (2000) and Kobayashi & Sari (2001) by assuming wind Lorentz 
factor distributions extending from r.^ ~ 10 to Tm ^ 10^. In our opinion the minimum wind 
Lorentz factor must be significantly higher than 10 (as otherwise the wind becomes optically 
thick) while the maximum Lorentz factor can not significantly exceed 10^ (due to the acceleration 
process limitations). 

We have shown that there is an upper limit to the observed energy of photons, Sp, at which 
the 7-ray flux peaks, Sp ^ e™^^ ^ O.BRqj^^ MeV (where Rq = lO^iioj cm), with very weak 
dependence on and on tv [see discussion preceding Eq. (^), and figures 2, 5 and 6]. Thus, the 
source dynamical time ~ Rq/c must satisfy ^ 1 ms in order to allow ~ 1 MeV. Sp ^ e™^^ 
is obtained for Tm ~ T^i, for which the radiative efficiency is largest, and Sp values in the range 
of 0.1 MeV to 1 MeV are obtained for wind parameters for which the radiative efficiency is high, 
^ 10% (see figure 3). 

High radiative efficiency and peak emission energy consistent with observations are therefore 
obtained for ~ Lmi- This does not necessarily imply that fine tuning of this model parameter 
is required. For F^ < Fm±, most efficient collisions occur at radii where the optical depth is high, 
leading to low efficiency, and hence low luminosity, bursts with peak emission energy ~ 1 keV (see 
figures 2 and 3), which would not have been detected by BATSE. For F^ significantly higher than 
Tm±, LF variance is small, leading to low efficiency, low luminosity bursts with peak emission 
energy ~ 10 keV, which may have been difficult to detect with BATSE. Natural consequences of 
the model considered here are therefore absence of bursts with peak emission energy significantly 
exceeding ~ 1 MeV, and existence of low luminosity bursts with low, 1 keV to ~ 10 keV, peak 
emission energy. The frequency of such bursts depends on the distribution of F^ in different 
winds. 

It should be pointed out in this context that while the deficit, among bursts detected by 
BATSE, of bursts with peak emission energy below 50 keV clearly reflects a real deflcit of such 
bursts, the deficit in detected bursts with peak emission energy > 0.5 MeV may be partly due to 
a selection effect (e.g. Llyod &: Petrosian 1999| ). BATSE data are consistent with equal number of 
bursts per logarithimic peak energy interval beyond Ep ^ 0.5 MeV (with luminosity comparable to 
that of lower peak energy bursts). Our model prediction, that bursts with ^ 1 MeV should be 
absent, should therefore be tested with future GRB detectors, which are able to better constrain 
the high energy end of the Sp distribution. 

A note should be made here regarding the low energy spectral slope. The upper limit on 
Tm is a consequence of the fact that for rj > r]^ shell acceleration saturates at F ~ r/,, < r/, as 
the wind becomes optically thin to Thomson scattering, and the internal energy left in the shell 
escapes as thermal radiation rather than being converted to kinetic energy. The requirement 
F ~ r/^, therefore implies that a significant fraction of the wind energy may escape as thermal 
radiation, leading to low energy spectral slopes steeper than those expected for pure synchrotron 
emission. This may account, at list partially, for observed steep low energy spectra. 
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Finally we note that the lower limit imposed on F^, ^ r±, is not derived from the 
requirement that the pair production optical depth for high energy, > 100 MeV, photons be 
smaller than unity. While this requirement leads to a similar constraint, ^ 10^, high energy 
photons have been detected in a small number of cases only. The constraint F^ > F^i is imposed 
in the present analysis by the requirement that the wind Thomson optical depth due to pairs 
produced by synchrotron photons be smaller than unity at the internal shocks stage. 
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Fig. 1. — Radiative efficiency versus peak emission energy for various Lorenz factor distributions. 
Different points within eacfi panel correspond to different choices of Tm, Tm and t^. Results are 
shown for 10 < Fm < Fm < 2500 and 10^"^ s < < 1 s, for two extreme assumptions on shell width 
expansion, maximal expansion A = max(i?0) -R/r^) and constant width A = Rq. We have assumed 
equal mass shells, Rq = 10^ cm, = lO^^ergs"^ and source redshift 2; = 1 for all calculations. 
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Fig. 2. — Contour plots of the peak emission energy, logiQ(ep[eV]), as a function of and the 
minimum Lorentz factor Tm of the ejected shehs. The value of Fm is 2500 and = 10^^ erg s~^, 
the LF distribution is bimodal and the shells have equal masses. The left panel refers to the case 
of expanding shell A = max(i?o, -R/F^) and the right panel to the case of time indipendent shell 
thickness A = Rq. 



-17- 




Fig. 3. — Contour plot of the radiative efficiency for tfie case described in Fig. 2. The radiative 
efficiency is similar for both A = max(i?o, R/T^) and A = Ro cases. 
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Fig. 4. — Same as right panel of Fig. 2 (constant pre-shock shell width), for electron and magnetic 
field energy fractions below equipartition. Left panel: Cg = 0.01, = 0.1; Right panel: Cg = 0.45, 
EB = 0.001. 
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Fig. 6.— 



Same as Fig. 2, for = 10^^ erg s ^. 



